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We apply novel jet techniques to investigate the spin and CP quantum numbers of a heavy 
resonance X, singly produced in pp — > X — > ZZ — > £ + £~jj at the LHC. We take into account all 
dominant background processes to show that this channel, which has been considered unobservable 
until now, can qualify under realistic conditions to supplement measurements of the purely leptonic 
decay channels X — > ZZ — ¥ 4£. We perform a detailed investigation of spin- and CP-sensitive angular 
observables on the fully-simulated final state for various spin and CP quantum numbers of the state 
X, tracing how potential sensitivity communicates through all the steps of a subjet analysis. This 
allows us to elaborate on the prospects and limitations of performing such measurements with the 
semihadronic final state. We find our analysis particularly sensitive to a CP-even or CP-odd scalar 
resonance, while, for tensorial and vectorial resonances, discriminative features are diminished in 
the boosted kinematical regime. 
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I. INTRODUCTION 

An experimental hint pointing towards the source of 
electroweak symmetry breaking (EWSB) remains miss- 
ing. From a theoretical perspective, unitarity constraints 
do force us to expect new physics manifesting at en- 
ergy scales < 1.2 TeV. For most of the realistic new 
physics models, this effectively means the observation of 
at least a single new resonant state in the weak boson 
scattering (sub)amplitudes, which serves to cure the bad 
high-energy behavior of massive longitudinal gauge bo- 
son scattering. An equal constraint does also follow from 
diboson production from massive quarks, qq — > W + W~ , 
which relates EWSB to the dynamics of fermion mass 
generation [1]. Revealing the unitarizing resonance's 
properties, such as mass, spin and CP quantum num- 
bers, is a primary goal of the CERN Large Hadron Col- 
lider (LHC) [2, 3]. Completing this task will provide 
indispensable information elucidating the realization of 
EWSB and will therefore help to pin down the correct 
mechanism of spontaneous symmetry breaking. 

At the LHC, the experimentally favorite channel 
to determine CV and spin of a new massive parti- 
cle X (mx > 300 GeV), coupling to weakly charged 
gauge bosons, is its decay to four charged leptons via 
X -> ZZ [4-8] . These "golden channels" [9] are charac- 
terized by extraordinarily clean signatures, making them 
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experimentally well-observable at, however, small rates 
due to the small leptonic Z branching ratios. Until 
now, X — > ZZ for the semihadronic decay channels 
ZZ — > £ + £~jj has been considered experimentally dis- 
favored. In part, this is due to the overwhelmingly large 
backgrounds from underlying event and quantum chro- 
modynamics (QCD), which exceed the expected signal by 
a few orders of magnitude. An analogous statement has 
been held true for associated standard model Higgs pro- 
duction, until, only very recently, new subjet techniques 
have proven capable of discriminating signal from back- 
ground for highly boosted Higgs kinematics [10-12]. Sub- 
sequently, these phase space regions have received lots 
of attention in phenomenological analysis, demonstrat- 
ing the capacity of subjet-related searches at the LHC in 
various channels and models [14]. H — > WW was one of 
the first channels subjet techniques were ever applied to 
[13]. The purpose of these analysis was predominantly to 
isolate a resonance peak from an invariant mass distribu- 
tion. In this paper we use a combination of different sub- 
jet methods for a shape-analysis of spin- and CP-sensitive 
observables in a semihadronic final state. CP properties 
of the standard model (SM)-like associated Higgs pro- 
duction have been investigated recently in Ref. [16] using 
^-tagging in HZ -> t+^bb. 

Subjet techniques have also proven successful in pp — > 
X — > ZZ — > £ + £~jj for standard model-like Higgs pro- 
duction (X = H with m„ > 350 GeV) in Ref. [15], 
yielding a 5<r discovery reach for an integrated lumi- 
nosity of 10 fb _1 , when running the LHC at a center- 
of-mass energy of s/s = 14 TeV. Equally important, 
the semihadronic decay channel exhibits a comparable 
statistical significance as the purely leptonic channels 
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pp —> X —> ZZ -» lllHi! , This can be of extreme im- 
portance if the LHC is not going to reach its center-of- 
mass design-energy. Hence, there is sufficient potential 
to revise semihadronic decays, not only to determine the 
resonance mass, but also its spin- and CV properties. 

The main goal of this paper is to investigate the at- 
tainable extent of sensitivity to the spin and CV quan- 
tum numbers of a resonance X in the channel pp — > X — > 
ZZ —> £ + £~jj, for the selection cuts, which allow to dis- 
criminate the signal from the background. To arrive at a 
reliable assessment, we take into account realistic simula- 
tions of both the signal and the dominating background 
processes. We fix the mass and the production modes 
of X, as well as its production cross section to be sim- 
ilar to the SM Higgs boson expectation*. On the one 
hand, this approach can be motivated by again referring 
to unitarity constraints: Curing the growth of both the 
VV — > VV and qq — > WW scattering amplitudes by a 
singly dominating additional resonance fixes the overall 
cross section to be of the order of the SM (see e.g. [17, 18] 
for non trivial examples). On the other hand, we would 
like to focus on an experimental situation, which favors 
the SM expectation, but leaving CV and spin properties 
as an open question. For this reason, we also do not in- 
clude additional dependencies of the cross section on the 
width of X. The width is, in principle, an additional, 
highly model-dependent parameter, which can be vastly 
different from the SM Higgs boson width (e.g. in models 
with EWSB by strong interactions [19, 20] or in so-called 
hidden- valley models [21]). Instead, we straightforwardly 
adopt the SM Higgs boson width, which then turns the 
resonance considered in this paper into a "Higgs look- 
alike" , to borrow the language of Ref . [8] . 

We organize this paper as follows: In Sec. II, we out- 
line the necessary technical details of our analysis. We 




FIG. 1: Spin- and CP-sensitive angles of Ref. [25] in pp — > 
X — > ZZ — > [i + fj,~jj. Details on the angles' definition and on 
the assignment of j a and jp are given in the text. An angle 
analogous to <3>i can be defined with respect to the leptonic 
decay plane. We refer to this angle as $. 



review the effective interactions, from which we com- 
pute the production and the decay of the resonance X 
with quantum numbers J CP = 0,1, 2 + . We also com- 
ment on the signal and background event generation and 
the chosen selection criteria, and we introduce the CV 
and spin-sensitive observables and their generalization to 
semihadronic final states. We discuss our numerical re- 
sults in Sec. Ill; Sec. IV closes with a summary and gives 
our conclusions. 



II. DETAILS OF THE ANALYSIS 

A. Spin- and CP-sensitive observables 

The spin and CV properties are examined through cor- 
relations in the angular distributions of the decay prod- 
ucts. A commonly used (sub)set of angles is given by 
the definitions of Cabibbo and Maksymowicz of Ref. [22] , 
which originate from similar studies of the kaon sys- 
tem (see, e.g., Refs. [4, 8, 23, 24] for their application 
to the X — > ZZ). In this paper we focus on the an- 
gles of Ref. [25] as sensitive observables, which also have 
been employed in the recent X — > Al investigation in 
Ref. [7]. We quickly recall their definition with the help 
of Fig. 1: Let p a , p^, and p± be the three-momenta of 
the (sub)jets j a and jp and the leptons in the laboratory 
frame, respectively. From these momenta, we compute 
the three-momenta of the hadronically and leptonically 
decaying Z bosons 

PZ h = Pa + Pf3 , PZ e =P++P , (la) 

as well as the lab-frame X three-momentum 



PX = Pa + PfS 



(lb) 



In addition, we denote the normalized unit vector along 
the beam axis measured in the X rest frame by e z , and 
the unit vector along the ZZ decay axis in the X rest 
frame by e z i. The angles of Fig. 1 are then defined as 
follows 



cos tih 



Pa ' PX 




, cos( 



P • Px 



P- Px 



(lc) 



z, 



(e z xe z ,)- (p_ x p + ) 



cos$ = 



, cos<fr— 

x v(P- x P+) 2 



(Pa X pp) ■ (p_ X P+ ) 

V(P* x Pp) 2 (P- x p+) 2 



X 

(Id) 

(le) 



"We normalize the cross section to SM Higgs production at the 
parton level. 



where the subscripts indicate the reference system, in 
which the angles are evaluated. More precisely, the he- 
licity angles 8h and 9g arc defined in their mother-Z's 
rest frame, and all other angles are defined in the rest 
frame of the particle X, where pz e = ~Pz h - It is also 
worth noting that the helicity angles correspond to the 
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so-called Collins-Soper angle of Ref. [26], evaluated for 
the respective Z boson. 

There is a small drawback when carrying over the defi- 
nitions of Eq. (1) from the purely leptonic decay channels 
to the considered semihadronic final state: When deal- 
ing with X — > £ + £~£' + £'~ , it is always possible to unam- 
biguously assign a preferential direction for the lepton 
pairs by tagging their charge^. This allows us to fix a 
convention for the helicity angles, as well as for the rel- 
ative orientation of the decay planes via a specific order 
of the three-momenta when defining the normal vectors 
in Eq. (1). Considering semihadronic X decays, we are 
stuck with a twofold ambiguity, which affects the angular 
distributions. Even worse, pr-ordered hard subjets, dug 
out from the "fat jet" during the subjet analysis (for de- 
tails see Sec. II C) can bias the distributions. Hence, we 
need to impose an ordering scheme which avoids these 
shortcomings. An efficiently working choice on the in- 
clusive parton level is provided by the imposing rapidity- 
ordering 

y(ja) < y{jp) , (2) 

which is reminiscent of the CP-sensitive <F,j observable 
in vector boson fusion [28]. This choice, however, docs 
not remove all ambiguities. The orientation of the de- 
cay planes [Eq. (le)] is not fixed by ordering the jets 
according to Eq. (2). The unresolved ambiguity results 
in averaging cos $ and cos(7r — $) over the event sample, 
leaving a decreased sensitivity in the angle $ <G [0, n]. We 
discuss this in more detail in Sec. III. 



B. Simulation of signal and background events 

We generate signal events pp — > X — > ZZ — > £ + £~jj 
with MadGraph/MadEvent [29], which we have 
slightly modified to fit the purpose of this work. In par- 
ticular these modifications include supplementing addi- 
tional Helas [30] routines and modifications of the Mad- 
GRAPH-generated code to include vertex structures and 
subprocesses that are investigated in this paper. We have 
validated our implementation against existing spin cor- 
relation results of Refs. [7, 32]. We choose the partonic 
production modes to be dependent on the quantum num- 
bers of the particle X: 

X = ± : gg -> X -> ZZ -> £ + £rjj , (3a) 
X = 1 ± : qq -> X -> ZZ -> £ + £rjj , (3b) 
X = 2+ : gg -> X -> ZZ -> t+T jj , (3c) 



t Charge tagging can be considered ideal for our purposes, given the 
additional uncertainties from parton showering; see Sec. III. For 
the region in pseudorapidity and transverse momentum that we 
consider, the mistagging probability of, e.g., muons is typically at 
the level of 0.5%, even for early data-taking scenarios [27]. 



where g denotes the gluon and q, j = (u, d, s, c) represents 
the light constituent quarks of the proton. 

The bottom quark contributions are negligibly small. 
While, in the light of the effective theory language of 
Ref. [7] , this specific choice can be considered as a general 
assumption of our analysis, the partonic subprocesses of 
Eq. (3) reflect the dominant production modes at the 
LHC. In particular, the production of an uncolored vector 
particle 1~ from two gluons via fermion loops is forbidden 
by Furry's theorem [33], while a direct ggZ' coupling is 
ruled out by Yang's theorem [34]. 

The effective operators that we include for the produc- 
tion and the decay of X do not exhaust all possibilities 
either (see again Ref. [7] for the complete set of allowed 
operators). Yet, we adopt a general enough set of op- 
erators to adequately highlight the features of objects 
X with different spins and CP quantum numbers in our 
comparative investigation in Sec. III. The effective ver- 
tex function, from which we derive the effective couplings 
of X to the SM Z bosons, that appear in the calculation 
of the matrix elements in Eq. (3), reads for the scalar 
case suppressing the color indices [35], 

£l» X =49^ + ^r ^vps Pi P2 ■ ( 4a ) 

For a vectorial X, the vertex function follows from the 
generalized Landau- Yang theorem [36] 

-c v 2 e^ pS (pi-p S 2 ), (4b) 
while for tensorial X 1 we include the vertex function [31] 

ZpvpS = C 1 (Pl,vP2, P g f ,S +Pl,pP2,p9vS 

+ Pi, P P2.8 9^v ~\ m x 9p, P 9p.s) (4c) 

to our comparison. 

From Eqs. (4), we can determine the (off-shell) decays 
X, X p , X p s — > Zp{pi)Z u (p2) by contracting with the final 
state Z bosons' effective polarization vectors £*Ap\) and 
£ ^(P2), which encode the Breit-Wigner propagator and 
the respective Z decay vertex. We include the spin and 
CV dependence of the X production from quarks via the 
effective Lagrangian in the vectorial scenario [30] 

L mX = ^y. {g v p L + g V R Vr) ^ ; (5a) 

where 

Vl,r= ^(ITTS) , (5b) 

project to left- and right-handed fermion chirality as 
usual. Defining g\ 2 — 9r± 9l> we can steer the vec- 
torial and axial couplings via g\ 2 • For the gluon- induced 
production of the scalar X case in Eq. (3), we compute 
the interaction vertices from 

£" X = -\ (.9? G^G^ X + g s 2 G^Gp V x) , (5c) 
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TABLE II: Cut flow comparison of the MadEvent signal event when processed either with Pythia 6.4 (referred to as P) or 
HerwigH — h (denoted by H) for the X states of Table I. Starting from the showered sample on calorimeter level (Raw), we apply 
the selection cuts (Cuts), the hadronic Z reconstruction requirements, the X mass reconstruction (mx), the S'/B-improving 
requirement on ARzz, and trimming and pruning (Tr+Pr). The selection criteria are described in detail in Sec. II C. The 
lower row gives the total cross section determined from the Herwig++ efficiencies after all cuts have been applied, cf. [15]. 
The different acceptance levels follow from the different angular correlations, see Sec. II C. 



Ref. [29] , where is the Hodge dual of the non- Abelian 
SU(3) field strength G^ v . For the production of the ten- 
sor particle X from gluons, we again assume the vertex 
function quoted in Eq. (4c). This choice corresponds to 
gravitonlike coupling, which, when taken to be universal, 
is already heavily constrained by Tevatron data (see e.g. 
[37] for recent D0 searches). The X = 2+, however, still 
represents a valid candidate for our spin and CV analysis 
as a state analogous to the composite A baryon. 

In the following we consider the five scenarios of Ta- 
ble I for our comparison in Sec. Ill for X mass and width 
choices m x = 400 GeV and T x = 27 GcV. The par- 
ton level Monte Carlo results for observables of Eq. (1) 
are plotted in Figs. 2-6 of Sec. III. From a purely phe- 
nomenological point of view, our strategy to normalize 
the parton level cross sections to the SM Higgs produc- 
tion at next-to-leading order (NLO, 50.84 fb after selec- 
tion cuts)'' 1 , effectively removes the dependence on the 
process-specific combinations of the parameters c^cf 
and gl'2, 9\, as well as the dependence on the initial state 
parton distribution functions on the considered spin- and 
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TABLE I: Definition of the scenarios considered for the com- 
parison in Sec. Ill for the X mass and width choices mx = 
400 GeV and T x = 27 GeV. 



■l-For the NLO Higgs production normalizations we use the codes 
of Refs. [38] and [39] for the gluon-fusion and weak-boson-fusion 
contributions, respectively. 



CP-sensitive angles. At the same time, the distinct angu- 
lar correlations will induce different signal efficiencies for 
the different particles X = J cv ', when the signal sample 
is confronted with the subjet analysis' selection cuts. In 
our approach, these naturally communicate to the final 
state after showering and hadronization. 

In principle, the s-channel signal adds coherently to 
the continuum ZZ production and their subsequent de- 
cay. We consider these Feynman graphs as part of the 
background and discard the resulting interference terms, 
which is admissible in the vicinity of the resonance. We 
have explicitly checked the effect of the interference on 
the angular distributions at the parton level for inclusive 
generator-level cuts and find excellent agreement for the 
invariant X mass window around the resonance, which is 
later applied as a selection cut in the subjet analysis. 

We further process the MADEvENT-generated signal 
events with HerwigH — h [40] for parton showering and 
hadronization. HerwigH — h includes spin correlations to 
the shower, hence minimizing unphysical contamination 
of the backgrounds' angular distribution by simulation- 
related shortcomings. We have also compared the results 
to Pythia 6.4 [41] to assess the systematic uncertainties 
and find reasonable agreement for the net efficiencies af- 
ter all analysis steps have been carried out (see Table II 
and the discussion of the next section) . 

We include the dominating background processes to 
our analysis. These arc Z+jcts, ZZ, tt and WZ pro- 
duction. For the main background Z + jets the NLO 
QCD cross section, requiring pr(jet) > 100 GeV, is 33.9 
pb [42]. The NLO QCD normalization of the tt produc- 
tion cross section is 875 pb [43], and for WZ production 
we find 43.4 pb [42]. The NLO QCD ZZ background is 
19.0 pb. For the simulation of the SM Higgs boson we 
take the NLO gluon-fusion and weak-boson-fusion pro- 
duction mechanism into account and normalize all sig- 
nals to the inclusive production cross section. We simu- 
late the backgrounds with MadEvent, again employing 
Pythia 6.4 and Herwig++ for showering, and normal- 
ize the distributions to the NLO QCD cross section ap- 
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FIG. 2: Cosine of the helicity angle 9h, Eq. (lc), calculated from the hadronically decaying Z at different steps of the analysis: 
inclusive Monte Carlo generation level (top, left), Monte Carlo generation level including selection cuts Eqs. (6)-(8) (top, right), 
after the full subjet analysis including Monte Carlo-truth information (bottom, left), and after the full analysis (bottom, right). 



proximation. For the weak-boson-fusion channels, the 
NLO QCD corrections are known to be at the percent 
level for a deep inelastic scattering type of factorization 
scale choices (see Ref. [39, 44]). For these processes, the 
electroweak modifications then become important since 
they turn out to be comparable in size [45]. As we per- 
form simulations within an effective electroweak theory 
approach, and for the smallness of the overall corrections, 
compared to additional uncertainties presently inherent 
to showering [46], we set a N /a ho = 1 for the weak- 
boson-fusion contributions. It turns out that the WZ 
background is negligibly small, and we therefore only fo- 
cus on the Z + jets, ZZ , and ti backgrounds in the dis- 
cussion of our results in Sec. III. 

From the fully-simulated Monte Carlo event, we recon- 
struct the detector calorimeter entries by grouping all fi- 
nal state particles into cells of size A77 x Atp = 0.1 x 0.1 in 
the pseudorapidity-azimuthal angle plane, to account for 
finite resolution of the calorimetry. The resulting cells' 
three-momenta are subsequently rescaled such to yield 



massless cell entries p^p^ = (see, e.g., Ref. [3]). For the 
rest of the analysis we discard cells with energy-entries 
below a calorimeter threshold of 0.5 GeV. 



C. Discriminating signal from background 

To separate the signal events of Eq. (3) from the back- 
ground, we perform a fat jet/subjet analysis. The tech- 
nical layout has been thoroughly discussed in the recent 
literature, e.g. in Refs. [10-12, 15]. We focus on the 
following on muon production. 

Prior to the jet analysis, we therefore require two iso- 
lated muons with 

Pt(m ± ) > 15 GeV and |??(m ± )| < 2.5 (6) 

in the final state, which reconstruct the invariant Z mass 
up to 10 GeV, 

[p(fi + )+p(fi-)] 2 = (m z ± 10 GeV) 2 . (7) 
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FIG. 3: Cosine of the helicity angle 9i, Eq. (lc), calculated from the hadronically decaying Z at different steps of the analysis: 
inclusive Monte Carlo generation level (top, left), Monte Carlo generation level including selection cuts Eqs. (6)-(8) (top, right), 
after the full subjet analysis including Monte Carlo-truth information (bottom, left), and after the full analysis (bottom, right). 



To call a muon isolated we require that i?T had < 0.1 Et ± 
within a cone of R = 0.3 around the muon. We ask for a 
"fat jet" with 

p T (fat jet) > 150 GeV and |y(fat jet) | < 2 , (8) 

defined via the inclusive Cambridge-Aachen algo- 
rithm [47] with resolution parameter R = 1.2. This par- 
ticular choice of R guarantees that we pick up the bulk 
of the Z decay jets since the boosted Z — ► jj decay prod- 
ucts' separation can be estimated to be ~ 2mz /pt(Z). 
Throughout, we invoke the algorithms and C++ classes 
provided by the Fast Jet framework [48]. 

The hadronic Z reconstruction is performed applying 
the strategy of Ref. [10]: For the hardest jet in the event 
we undo the last stage of clustering, leaving two subjets, 
which we order with respect to their invariant masses 
rrij 1 > nij 2 . Provided a significant mass drop for a not 
too asymmetric splitting, 

m h < 0.67 mj , AR 2 J1J2 min(^ , ^ ) > 0.09 m) , 

(9) 



where AR denotes the distance in the azimuthal angle- 
pseudorapidity plane, we consider the jet j to be in the 
neighborhood of the resonance and terminate the declus- 
tering. Otherwise we redefine j to be equal to ji and con- 
tinue the algorithm until the mass-drop condition is met. 
In case this does not happen for the considered event, 
and we discard the event entirely. If the mass-drop con- 
dition is met, we proceed with filtering of the fat jet [10]; 
i.e. the constituents of the two subjets which survive the 
mass drop condition are recombincd with higher resolu- 
tion 

fl mt =mm(o.3,^p) , (10) 

and the three hardest filtered subjets are again required 
to reproduce the Z mass within mz ±10 GeV. 

We subsequently reconstruct the Higgs mass from the 
excess in the m 2 x = (pz h + Pz e ) 2 distribution, i.e. we 
imagine a situation where the X mass peak has already 
been established experimentally. This allows us to avoid 
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FIG. 4: Cosine of the angle 6*, Eq. (Id), calculated from the hadronically decaying Z at different steps of the analysis: inclusive 
Monte Carlo generation level (top, left), Monte Carlo generation level including selection cuts Eqs. (6)-(8) (top, right), after 
the full subjet analysis including Monte Carlo-truth information (bottom, left), and after the full analysis (bottom, right). 



dealing with the sophisticated details of experimental 
strategies, which aim to single out the resonance peak 
from underlying event, pile-up and background distri- 
butions by typically including combinations of various 
statistical methods . A thorough discussion would be be- 
yond the scope of this work. For a considered X mass 
of 400 GeV we include events characterized by recon- 
structed invariant masses 

[p(/i+) +p(jT) + P (j a )+P(jp)] 2 = (400 ± 50 GeV) 2 . 

(11) 

Further signal-over-background (S/B) improvements can 
be achieved by requiring 

AR ZZ < 3.2 

and by trimming and pruning [49] of the hadronic Z event 
candidates on the massless cell level of the event [12], as 
described in detail in Ref. [15]. 



III. RESULTS AND DISCUSSION 

In Figs. 2-7 we show the angles of Eq. (1) after various 
steps of the analysis have been carried out. We also give 
a comparison of the full hadron-level result and Monte 
Carlo truth; i.e. we take into account the shower's parti- 
cle information. These plots are the main results of this 
paper. Comparing the two shower and hadronization ap- 
proaches of Pythia 6.4 and Herwig++ for the process 
efficiencies in Table II, we find substantial discrepancies 
at intermediate steps of our analysis. After the entire 
analysis has been carried out this translates into a sys- 
tematic uncertainty of ~ 30% of the total cross sections. 
This is not a too large disagreement as both programs 
rely on distinct philosophies and approaches, which typi- 
cally result in sizable deviations when compared for iden- 
tical Monte Carlo input. The plots in Figs. 2-7 show 
distributions obtained with HerwigH — h 

We now turn to the discussion of the angular correla- 
tions. It is immediately clear that the chosen selection 
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full reconstruction, Herwig + + 



FIG. 5: Angle $, Eq. (Id), calculated from the hadronically decaying Z at different steps of the analysis: inclusive Monte Carlo 
generation level (top, left), Monte Carlo generation level including selection cuts Eqs. (6)-(8) (top, right), after the full subjet 
analysis including Monte Carlo-truth information (bottom, left), and after the full analysis (bottom, right). 



criteria, Eqs. (6)-(8), do heavily affect the sensitive an- 
gular distributions of Eqs. (1). Retaining a signal-over- 
background ratio of approximately 0.5, however, does not 
allow us to relax the pt cut on the fat jet. This cut 
turns out to be lethal to some of the angular distribu- 
tions. Referring, e.g., to cos 9*, plotted in Fig. 4, we find 
that our fat jet criteria, Eq. (8), force the distribution 
into reflecting extremely hard and central decay prod- 
ucts. This removes essentially all discriminating features 
from the differential distribution dcr/dcos#*, that show 
up for | cos 0*| > 0.5 at the (inclusive) Monte Carlo event 
generation level. This is also reflected in the distinct 
acceptance level of the different J samples, shown in 
Table IF Note that, throughout, the fully hadronic dis- 
tributions are in very good agreement with the Monte 
Carlo-truth level. 

Most of the sensitivity found in the observable $ for 
the signal sample can be carried over to the hadron level. 
Yet, the angular pattern is known to be sensitive to the 
X's mass scale, tending to decorrelate for larger X masses 



(see e.g. [7]). 

As already pointed out in Sec. II, the ambiguity in 
cos $ smears out the angular correlations quite a lot in 
Fig. 5. This comes not as too large limitation of the 
angle's sensitivity for a CP-odd scalar particle X. For 
X = 0~ , the distribution peaks at <E> = ir/2 and is also 
rather symmetrical with respect to tt/2. This leaves us af- 
ter the subjet analysis with the helicity angles of Eq. (lc) 
and $ as three sensitive angles out of five not taking into 
account the background distribution. 

Crucial to obtaining angular correlations after all, is 
the analysis' capability to reconstruct both of the Z rest 
frames (and from them the X rest frame). This is al- 
ready clear from the angles' definition in Eq. (1), and, 
again, this is not an experimental problem considering 
the purely leptonic channels. For the angles $ and 0* 
decorrelate (with the exception of 0~) due to the selec- 
tion criteria, a bad rest frame reconstruction would not 
be visible in these observables immediately. This is very 
different if we turn to the helicity angles. Quite obviously, 
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FIG. 6: Angle <I>, Eq. (le), calculated from the hadronically decaying Z at different steps of the analysis: inclusive Monte Carlo 
generation level (top, left), Monte Carlo generation level including selection cuts Eqs. (6)-(8) (top, right), after the full subjet 
analysis including Monte Carlo-truth information (bottom, left), and after the full analysis (bottom, right). 



given a good hadronically decaying Z rest frame recon- 
struction, we can apply the identical leptonic helicity an- 
gle as invoked for the measurement in X —> ZZ —> A£; we 
have referred to this angle as 0£ , previously. The only dif- 
ference compared to the purely leptonic analysis is that 
we consult a partly hadronic system to construct the ref- 
erence system, in which the leptonic helicity angle Qg is 
defined. 

Indeed, the subjet analysis described in Sec. II C is ca- 
pable of giving a very good reconstruction of the hadron- 
ically decaying Z boson rest frame, while sufficiently re- 
ducing the backgrounds. This allows to carry over most 
of the central sensitivity of the angular distributions in 
Fig. 3 to the fully simulated final state. However, the 
hadronically-defined helicity angle, displayed in Fig. 2, 
also suffers badly from the subjet analysis. Note that 
the bulk of the modifications of cos Oh do not arise from 
our restrictive selection criterion Eq. (8), but from sym- 
metry requirements among the subjets in the mass-drop 



procedure. Thus, the subjets which provide a significant 
mass drop with respect to Eq. (9) are biased towards 
9 h -90°. 

A remaining key question that needs to be addressed is 
whether the potentially sensitive angles 6*£,<I>, and $ ex- 
hibit visible spin- and CP-dependent deviations when the 
background distribution is taken into account. We show 
these angles including the backgrounds in Fig. 7. The 
backgrounds' $ distribution largely mimics the 1 + shape 
under the subjet analysis' conditions, so we cannot claim 
sensitivity unless the backgrounds distribution is very 
well known. This also accounts for the $ distribution in 
a milder form. While here the background is flat to good 
approximation, S/B (see Table II) limits the sensitivity 
to the shape deviations, which are ameliorated due to the 
different signal efficiencies. However, the distribution re- 
mains sensitive to X = CP shape, pp — > X — > ji + fj,~ jj 
remains sensitive to the CP quantum number of a scalar 
particle X in the cos 9g distribution, which is opposite in 



10 




0.5 1.0 1.5 2.0 2.5 

full reconstruction, Herwig + + 




0.5 1.0 1.5 2.0 2.5 3.0 

full reconstruction, Herwig + + 




-0.5 0.0 0.5 

cosS,, full reconstruction, Herwig + + 



FIG. 7: The spin- and CP-sensitive angles $ (top), 
backgrounds, simulated with Herwig+ + . 



$ (bottom, left) and cos9e (bottom, right) including the shape of the 



shape compared to the background distribution. 

We have only considered an X mass mx = 400 GeV, a 
choice which is quite close to the lower limit of the mass 
range, where the boosted analysis is applicable. Some re- 
marks concerning our analysis for different X masses and 
widths are due. The boost requirements and the centrally 
required selection cuts do affect the angular distributions 
in a X mass-independent manner. The remaining an- 
gles are then qualitatively determined by the goodness 
of reconstruction, which becomes increasingly better for 
heavier X masses, keeping the width fixed. 

In case of the SM Higgs boson, the width is propor- 
tional to m^ H due to the enhanced branching of the Higgs 
to longitudinally-polarized Zs. With the resonance be- 
coming width-dominated, our mass reconstruction still 
remains sufficiently effective; S/B, however, increasingly 
worsens. For these mass ranges, the analysis is sensi- 
tive to the experimental methods that recover the reso- 
nance excess. Additionally, from a theoretical perspec- 
tive, there are various models known in the literature 
where a heavy resonance becomes utterly narrow or ex- 



ceedingly broad (see, e.g. [18, 50] for discussions of the 
resulting phenomenology). The former yields, depend- 
ing on the (non-SM) production cross section, a bet- 
ter mass reconstruction, while the latter case is again 
strongly limited by S/B; cf. Fig. 7. For any of these 
EWSB realizations, our methods should be modified ac- 
cordingly, taking into account all realistic experimental 
algorithms, techniques, and uncertainties as well as all 
mo del- dependent parameters. 



IV. SUMMARY AND CONCLUSIONS 

The discovery of a singly produced new state at the 
LHC, and methods to determine its additional quantum 
numbers, remains an active field of particle physics phe- 
nomenology research. The availability of tools to simu- 
late events at a crucial level of realism and precision has 
put the phenomenology community into a position that 
allows a more transparent view onto the complicated par- 
ticle dynamics at high-energy hadron colliders than ever. 
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In this paper we have explored the performance of 
new jet techniques when applied to the analysis of spin- 
and CP-sensitive distributions of a newly discovered res- 
onance, which resembles the SM in the overall rate in 
pp — )■ We have performed a detailed investiga- 

tion of the angular correlations and have worked out the 
approach-specific limitations, resulting from the boosted 
and central kinematical configurations. It is self-evident 
that a QCD-dominated final state cannot compete with a 
lcptonic final state in terms of signal purity, higher order 
and shower uncertainties, per se. These uncertainties are 
inherent to any current discussion related to jet physics. 
Nonetheless, we have shown that potential "no-go the- 
orems" following from huge underlying event and QCD 
background rates for pp — > X — > ZZ — > l + £~jj can be 
sufficiently ameliorated to yield an overall sensitivity to 
the CP property of a singly produced scalar resonance. 
Straightforwardly applying the described analysis strat- 
egy to vectorial and tensorial resonances does not yield 
reliable shape deviations when the backgrounds' distribu- 
tion is taken into account. Given that the cross section 
of the semihadronic decay channel is approximately 10 
times larger compared to X — > At, the performed subjet 
analysis qualifies to at least supplement measurements of 
the purely leptonic decay channels. 

A question we have not addressed in this paper is the 
potential application of the presented strategy to sig- 
natures, which do not resemble the SM at all. Elec- 
troweak symmetry breaking by strong interactions is 
likely to yield a large rate of longitudinally polarized elec- 
troweak bosons due to modified XZZ and Xqq couplings 



[19, 20]. Measuring the fraction of longitudinal polariza- 
tions, which can be inferred from the Z's decay prod- 
ucts' angular correlation as proposed recently in Rcf. [51], 
should benefit from the methods we have investigated in 
this paper. This is in particular true for new compos- 
ite operators, such as a modification of the Higgs kinetic 
term [20], inducing asymmetric angular decay distribu- 
tions of the leptons. In addition, our analysis is also 
applicable to the investigation of isovectorial resonances 
(see, e.g., Ref. [50]) in pp — > WZ, with the W decaying 
to hadrons, and Z — > n + n~ ■ We leave a more thorough 
investigation of these directions to future work. 
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